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Q\ • The Sudden Approximation is applied to invert structural data on randomly corrugated surfaces 

0^ ' from inert atom scattering intensities. Several expressions relating experimental observables to sur- 

face statistical features are derived. The results suggest that atom (and in particular He) scattering 
(— i ■ can be used profitably to study hitherto unexplored forms of complex surface disorder. 
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\Q ; I. INTRODUCTION 

Structurally disordered surfaces have been a subject of great interest for some time now. Of special interest are 
^ ■ epitaxially grown films, liquid surfaces, and amorphous surfaces. In epitaxial growth for example, metal or semicon- 
ductor atoms are adsorbed on a surface under thermal conditions, to form two- and three-dimensional structures on 
top of it. The physical and chemical properties are determined by the final form of these structures. These may be 
of dramatic importance, e.g, in the production of electronic devices. One of the most exciting aspects of epitaxial 
growth kinetics, is that it prepares disordered structures in the intermediate stages. The disorder manifests itself in 
the formation of various types of clusters or diffusion limited aggregates on top of the surface. These structures may 
be monolayers (usually at high temperatures, when the diffusivity is large, or at coverages significantly below a mono- 
layer), in which case the disorder is two-dimensional, or they may be composed of several layers, giving rise to disorder 
in three dimensions. Epitaxially grown structures of this type offer an exceptional opportunity for both experimental 
and theoretical study of disorder. No satisfactory and comprehensive theory of the epitaxial growth process is as of 
yet available, much due to the absence of reliable interaction potentials for the system. The situation with respect 
to liquid and amorphous surfaces is similar: very little is known at this point about their structure. Progress at 
this stage thus hinges critically on data available from experiments. An important experimental technique is thermal 
atom scattering, and in particular He scattering The main advantage offered by He scattering is complete 

surface sensitivity, as He does not penetrate into the bulk, unlike other scattering techniques such as neutron or X-ray 
scattering, or low energy electron diffraction (LEED). Another important advantage is that He scattering is highly 
non-intrusive, due to the inertness and low mass of the He atoms. The latter also means that He scattering is really 
a diffraction experiment at the typical meV energy scale at which most experiments are performed, with sensitivity 
to atomic-scale features. The interpretation of He scattering experiments is, however, rather involved due to the 
| complicated interaction between the He beam and the surface. As in all other scattering problems, this interpretation 
issue is in fact one of inversion of the He/surface potential. The inversion problem, however, is intrinsically ill-posed, 
since one can only measure intensities, not phases. There is a certain redundancy in the intensity information which 
can be exploited to obtain relative phase shifts, but never the absolute phases. Thus the inversion can never be fully 
performed, although partial information can be obtained, or useful approximations can be made concerning the shape 
of the potential, which yield an analytically closed-form solution . 

In this work I will concentrate on the problem of trying to relate surface structural features to the scattering 
O ■ intensities, rather than the potential in general. This involves a non-trivial step of connecting the potential with such 
features. One possibility is the suggestion by Norskov and coworkers |7) on the basis of effective-medium theory, in 
which the leading repulsive term in the He scattering potential from any electronic system is taken to be proportional 
to the local unperturbed electron density no(r) of the host at the He position r: 
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V{r) = cm (r). (1) 

In this way the potential is rather simply related to structural features such as the local electronic corrugation. 
However, while useful for describing the interaction close to the surface, where repulsion dominates, this formula is 
inapplicable at large distances where the He/surface interaction is dominated by the long-range attractive forces. In 
addition, it still requires knowledge of the local electron density, often a highly nontrivial task, especially in presence 
of defects or impurities. Thus in practice one often resorts to the use of specific functional forms for the potential, 
such as the Morse or Lennard- Jones potentials, and fits the parameters to the experimental data. Unfortunately, the 
connection to surface geometrical features is then less transparent. 
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So far almost all of the work done on inverting He/surface potentials has been for He scattering from ordered, 
crystalline surfaces. Again, a full inversion is impossible, but useful results have been obtained by assuming specific 
forms for the interaction potential. These include semiconductor surfaces: GaAs(llO), Si, InSb(llO) |^|, and some 
transition metal surfaces (Ni, Cu, Ag, Au) 0). More recent efforts have concentrated on potentials for scattering from 
various adsorbed monolayers on metals, e.g., the c(2*2) phase of oxygen on Ni(OOl) [[To), hydrogen-plated Pt [jll], and 
Xe/Cu(110) ]l2] ]. We recently fitted cross-section data to obtain a potential for a disordered Ag/Pt(lll) system Jl3| . 
A very detailed recent review of potentials of physical adsorption is the work of Vidali et al. E3] , which tabulates 
parameters of interest as deduced from analysis of experimental data and calculations of over 250 gas-surface systems, 
including He. 

Formal inversion methods have also occupied the attention of various researchers. The first such work was presented 
by Gerber and Yinnon Q , who showed that atom-surface interaction potentials can be recovered from the diffraction 
peak intensities measured in beam scattering experiments by a direct, simple inversion method, using the Sudden 
approximation (SA). The SA is a highly successful and useful theoretical method in the He scattering field, and has 
been reviewed by Gerber The method of Ref. was applied to simulated Ne scattering from W(110). This was 
followed by the first inversion of real atom/surface scattering data: the He/MgO(100) system Schlup [|6| showed 
how the surface profile function can be inverted in the Eikonal approximation, and applied his method to simulated 
data. In Ref. H we studied a related problem: the inversion of an ad-atom profile function in the SA. Rabitz and 



coworkers developed a method based on functional sensitivity analysis and Tikhonov regularization |17|, and applied 
this to the inversion of the He potential for scattering from a Xe monolayer on the (0001) face of graphite |p|,^8|, which 
was also attempted earlier by employing close-coupling calculations Jl9[ |. Finally, we recently inverted the structure 
of a low temperature disordered overlayer of Ag on Pt(lll) po| . 

It should be noted that the situation with regards to potentials for He scattering from disordered surfaces is very 
much inferior to the case of ordered surfaces described above. Essentially no such reliable potentials exist, and 
the subject is at this point at a most preliminary stage. The common approach it to assume that the potential 
can be represented as a surface term plus a sum of pairwise additive terms representing the interaction of the He 
with each of the surface adatoms or vacancies. In this paper I will take a different approach, and will show how 
one can derive statistical information about randomly corrugated surfaces from He scattering measurements. This 
information comes in the form of correlation functions, not potential parameters. The resulting expressions relate 
useful statistical parameters characterizing the surface disorder to experimental observables which are not hard to 
obtain in practice, such as incidence energy dependence of the specular intensity. Unfortunately, at the time of 
writing experiments are unavailable for comparison with the results obtained here. The developments will therefore 
be primarily methodological, in anticipation of experimental data. It is hoped that the results obtained here will 
motivate He scattering experiments on disordered solid and liquid surfaces. As demonstrated by this work and others 
before it, He scattering can provide a wealth of information on disordered surface structure and d ynam ics. 

The structure of the paper is as follows. Sec. || provides a brief introduction into the SA. Sec. 



the paper and derives the inversion expressions. Concluding remarks are brought in Sec. IV 



III 



is the heart of 



II. BRIEF REVIEW OF THE SUDDEN APPROXIMATION 

Consider a He atom with mass /x incident upon a surface with wavevector k= (K, k z ). KK. — h(k x , k y ) and KK' are 
respectively the intial and final momentum components parallel to the surface, and hk z is the incident momentum 
normal to the surface. The position of the He atom is r— (R, z), where R= (x,y) is the lateral position. The SA is 
valid when the collisional momentum transfer q = K' — K in the direction parallel to the surface is much smaller than 
the momentum transfer normal to the surface: 2k z 3> |q|. This condition is satisfied at relatively high incidence angle 
and energy E = (Kk) 2 / (2m) , and moderate surface corrugations. When it holds, one can approximately consider the 
scattering along z as occurring at fixed R. Then if ?p is the He wavefunction, it satisfies a Schrodinger equation where 
the dependence on R is adiabatic: 



<Pr(z) = #r(z). (2) 



Here Vr (z) is the He-surface interaction potential and no inelastic channels are included, so that the total energy 
e is conserved . This means that each surface point R gives rise to an elastic real phase shift ?y(R), which can be 
evaluated in the WKB approximation from Eq.(||) as: 
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(3) 



where £(R) is the classical turning point pertaining to the integrand in Eq.(||). The phase shift in turn yields the 
S-matrix as: <S(R) = exp[2iry(R)]. The R coordinate is conserved in this picture so the S-matrix is diagonal in the 
coordinate representation: 

(R'|5|R) = e 24 " (R) <5(R' -R). 

Experimentally one measures probabilities |<S(K' — > K)| 2 for K — * K' transitions. To obtain these (R|K) = exp(iK- 
R)/^/A (where A is the area of the surface) can be used, to find: 



(K'|5|K) = J dR' dR(K'|R')(R'|5|R)(R|K) =jj dRe-' iq - R e 2l??(R) . 



(4) 



This is the well-known expression for the SA scattering amplitude plj . The SA has been tested extensively by 
comparison to exact coupled-channel calculations on Ne/W(110) and He/LiF(001) p2| , as well as by comparison to 
exact time-dependent propagation methods in the case of scattering from defects p3|7 It is particularly noteworthy 
that the SA yielding an inverted potential of remarkable accuracy Q for simulated Ne/W(110) data. The most 
important deficiency of the SA (apart from being limited to high energies) is its inability to describe double collision 
events. This is because, as mentioned above, the R coordinate is conserved in the SA, i.e., each trajectory takes place 
at constant R. Clearly, no double collisions can occur under such conditions. However, double collisions may take 
place, e.g., when an incident atom is scattered off a defect onto the surface, or in the opposite order. This issue was 
resolved recently by combining the SA with the Born approximation p4fl. 



III. INVERSION OF STRUCTURE OF RANDOMLY CORRUGATED SURFACES FROM He 

SCATTERING INTENSITIES 

A. Atom Scattering from a Randomly Corrugated Square Well Potential 

Recall the role of £(R) in the SA amplitude [Eq.(^)]: it is the position of the classical turning points, which can 
alternatively be viewed as the surface corrugation function. In the case of a disordered surface £(R) can be a random 
function of R and in order to obtain observable quantities one must average over an appropriate ensemble which 
characterizes the physical and statistical properties of the surface of interest. Thus the scattering probability, as a 
function of momentum transfer Q, is given by 

p(Q)H|S(Q)| 2 ), (5) 

where (• • •) indicates an average over the ensemble of which £(R) is a typical sample. In certain cases one will be 
justified in assuming that translational invariance has been established after averaging over the corrugation ensemble. 
For example, this will generally be the case for liquid surfaces, for solids when the structural disorder on the surface is 
due to radiation damage, and epitaxially grown defects on a surface. This means that £(R) is a stationary stochastic 
process, and hence that (/[£(R)]) is independent of R, and that (g[£(R), £(R')]) depends only on R — R'. The 
assumption of stationarity will, however, not hold for the case in which the structural disorder is caused by an 
incomplete adsorbed overlayer on a periodic substrate, and I will present a different treatment for such cases in a later 
section. Combining Eqs.(Q),(^) one obtains in the case of a translationally invariant potential 

P( Q ) = JL y dRe iQ-R( e HvOV-n(n)]^ (6) 

From the above expression it can be seen that within the SA, the angular scattering intensity is the Fourier transform 
of a function of the random variable 

E(R) = 2[r ] (0)-r l (R)}. (7) 

To proceed, one must first establish the connection between the phase-shift and properties of the surface. Assume 
now that the interaction of the He with the entire surface can be expressed in the form of a square well potential of 
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depth e. Square well potentials have been rather successful in predicting properties of interacting gases; the accuracy 
obtained in fitting second and third virial coefficients with a square well potential compares to that of a Lennard- Jones 
6-1 2 potential ]25}] . They have also been studied by others in the context of inversion problems [ p6|]27| , e.g., in neutron 



reflectometry from magnetic films 28 1. In our case the potential assumes the form: 

2m 
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_ 2 V(R,z)=^ -e 




z < £(R) 

£(R) < z < £(R) + A£(R) = C(R) 
z > C(R) 



Typically, £(R) is a relatively strongly corrugated function as it expresses the interaction of the He atom with the 
core electrons responsible for the steep repulsive part of the potential. On the other hand, £(R) m ay be very smooth, 
reflecting the loss of detail at the long distances at which the attractive part of the He-surface potential becomes 
important. According to Eq.(||) the phase shift is then given by 

7? (R) = k z A£(R) ((1 + S)i -l)-k z £(R), (8) 

where S = e/k z is a small parameter in the high-energy SA. Eq.(||) is known as "the Beeby effect" [M, i.e., in the 
presence of a well the wave number has to be replaced by an effective wave number which is due to the acceleration 
of the particle by the attractive forces. In the case of scattering by a pure hard wall [A£(R) = 0], or a square well 
with the same shape as the hard wall [A£ (R) = z = constant] , the scattering intensity is 

P(Q) = I J dRe iQ-R< e 2ifc*K(0)-e(R)]). (9) 

In this case the intensity is essentially the Fourier transform of the characteristic function exp[ik z Z(R)] of the relative 
surface corrugation: 

Z(R) = 2fc(0) - £(R)], 

which makes its interpretation particularly clear. Knowledge of the probability density /^(z;R) fully determines 
the corrugation of the surface in question, in that fz(z;TV)dz is the probability of finding a corrugation of height 
between z and z + dz at R, with respect to that at the origin. It is a well known fact in probability theory |3C| ] 
that for any (piecewise-) continuous distribution function fx{x) a unique and explicit inversion of fx( x ) exists from 
the characteristic function. Hence it is clear at the outset that He scattering can, at least in principle, provide very 
useful information about the statistics of a disordered system. In the following sections I will show how the statistical 
information pertaining to the disordered surface can be extracted from the He scattering intensity. 



B. Extraction of the Surface Corrugation Probability Density from the Angular Scattering Intensity for a 

Hard Wall Potential 

In this section I specialize to the hard wall potential. This model has been in use for many years in He scattering 
theory J3l]-|33| and is related to the Eikonal approximation in optics (34). In Ref. J32| the origin of the hard- wall 
was discussed and it was concluded that it is due to the the local density of metal electron states in the selvedge. 
This leads to the Esbjerg-Norskov theory, Eq. ([[]). The great merits of the hard wall approximation are that (1) it 
is analytically tractable and (2) it provides a direct geometrical interpretation of the surface corrugation. However, 
it is clearly oversimplified and leaves out many interesting features of the He-surface interaction. More sophisticated 
approximation schemes have therefore been introduced by various researchers. The distorted-wave Born approximation 
was used early on pq] to treat soft potentials |Q and low-energy scattering J37] . These studies have shown significant 
deviations from the predictions of the hard- wall model, e.g., for a corrugated Morse potential for He scattering by a 
Cu(llO) surface p^j. Nevertheless, the hard wall approximation has physical merit in the He high-energy limit, which 
is assumed here in connection with the SA. 

An inverse Fourier transform of the scattering intensity P(Q) [Eq.(^)] yields the characteristic function of the 
surface corrugation: 

(e i^(R) ) = ^l [p(Q);R]) (1Q) 
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where a notation for the Fourier transform was introduced: 

T[f{x)-y] = ^= r dxf(x)e ix v. 

In the case of normal incidence -P(Q) is a symmetric function, so that its Fourier transform is real. The probability 
density fz(z; R) can now be found if one recalls that: 

(e ik ° z ^) = J dz e lk > z f z (z; R). (11) 

Formally, therefore, /^(z;R) is fully determined by the scattering intensity, through another Fourier transform: 

f z (z;R) = J dk z e- ik **F^[P(Q);B). (12) 

The last equation shows that He scattering may in principle be used to fully invert the probability density of the 
corrugation function of a randomly corrugated surface. However, this requires dense sampling of the scattering 
intensities over a broad range of incidence wave-numbers, a task which may be difficult to accomplish in practice. 
Additional difficulties may arise due to Gibbs phenomenon and noise. A more realistic approach is to determine 
moments of the random variable describing the surface corrugation, an idea which dates back to the Backus and 
Gilbert work in geophysics J|^] , and has been pursued by others as well |2^,^] . To show how this may be accomplished 
here one may expand the characteristic function in terms of its moments. We then obtain from Eqs.(10),(|i"T|): 



E ^(W) = ^ X [P(Q);R] = 9(k z ;R). (13) 

The moments of Z(R) may now be found by differentiation of the experimentally available function g(k z ; R): 

(z(Rn = _LM. (14) 



Thus in contrast to the rather involved task of extraction of the full probability density via Eq.(12), requiring an 
experiment to be performed over the entire energy range, extraction of the moments merely involves an extrapolation 
of the data to the low energy limit k z = 0. Progressively higher moments, however, require higher derivatives and 
hence an increasingly dense sampling in k z to reduce noise due to finite differences. As a check, one may estimate 
the first moment (Z(R)) = 2(£(0) — £(R)), which must clearly vanish. In practice, perhaps the most interesting piece 
of information is the second moment. Assuming that all higher moments vanish is equivalent to the assumption that 



Z(H) is a Gaussian random variable. From Eqs.(|l3|),(14) one has: 



<Z(R)2> = Ki p mn]] k= . (15) 

One may also choose to focus on the corrugation function £(R) itself, instead of on the relative corrugation Z(R) = 
2[£(0) — £(R)]. Without loss of generality one can take (£(0)) = (£(R)) =0. A cumulant expansion can then be used 
to evaluate the averages in Eq. ([)]). Truncating this expansion at the second cumulant, or equivalently, assuming that 
£(R) is a Gaussian random variable, one obtains: 

p(Q) = _L I rfRel QR e (2 ff fe 2 ) 2 [C(R)-l] ; (16) 



where the variance and correlation function are respectively 

a 2 = (£(R) 2 ), C(R) = -L(£(0)£(R)). (17) 



Here, by definition, 0< |C(R)| < 1 and C(R) — > as R — > oo. Fourier transforming the scattering intensity now yields 
an expression which can be used to conveniently fit a 2 and C (R) as a function of k z : 

A 



e 



(2cr k z 



[C(R)-1] = il^-l[p ( Q );R] . (18) 
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A set of equations [notably Eqs.(Q),( 15 ) , (Jlq)] has thus been derived which can be applied to extract useful statistical 
information on the surface corrugation, within a hard wall model, from the angular intensity distribution. In the next 
section I will consider what information can be derived within the more realistic square well model. 
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C. Correlation Functions in the Square Well Model 



It is not possible to analytically solve for the probability density of the effective surface corrugation or its moments 



in the square well model, as was done in Sec. Ill B in the hard wall case. It is, however, possible to obtain some 
interesting information by performing a cumulant expansion, as I now proceed to show. An inverse Fourier transform 
of the scattering intensity in the general case of a translationally invariant potential [Eq.(^|)] yields [using Eq. (]?])]: 

(e iK(R) } = i.^-l[p(Q);R]. 

Introducing the notation 

AZ(R) = 2[A£(R) - A£(0)]. 
one obtains with the help of Eq.(|h, after some algebra: 



E(R) 



Z(R) - -8AZ(R) 



This expression is exact to first order in the small parameter S. Without loss of generality again set (Z(R)) =0, and 
also (AZ(R)) =0, so that the first moment of E(R) vanishes. For the second moment one obtains, again to first order 
in 5: 

(E(R) 2 ) - (E(R)) 2 » k 2 [(Z(R) 2 ) - 5{Z(R) AZ(R))] . 

In analogy to Eq. ([l7]) define a variance and correlation function between the hard wall corrugation and the deviation 
from it due to the attractive well: 

a\ = (£(R) A£(R)), Cx(R) = -i(£(R) A£(0)), 

a i 

so that 

(Z(R) AZ(R)} = 8a 2 [l - Ci(R)]. 
Collecting the results we obtain finally to second order in the cumulant expansion: 

e (2<rfc i ) 2 [C(R)-l]-(2<T 1 e) 2 [C 1 (R)-l] = A JT~l [P(Q); R] . 

This equation [compare to Eq.(|l8|)] may be used to fit the experimental data to obtain the correlation functions, 
variances, and the well depth of the attractive part of the potential. This is much information to fit to a single 
experimental function. However, since the contribution due to the attractive well is expected to be small, one may at 
first neglect this contribution. This is tantamount to assuming a hard wall interaction, by which one may obtain an 
estimate of the hard wall quantities a and C(R). Having computed these, one may correct the fit by inclusion of ci, 
Ci(R), and e, and proceed iteratively. In the next section I will consider how statistical information can be derived 
from the experimentally straightforward measurement of the specular intensity. 



D. Extraction of the Correlation Length from Specular Peak Measurements in the Hard Wall Model 

From an experimental point of view, there is a significant advantage to working with specular intensities, as 
their measurement involves a relatively minor effort compared to that required to obtain the full angular intensity 
distribution. There is also a theoretical advantage, in that energy transfer to the surface is significantly reduced in 
specular collisions, so that a phononless treatment is more justified. In the hard wall approximation, by addition and 
subtraction of the same term, Eq. (j^) for the angular distribution can be written as 



P(Q) 
1 

A* 



dRe ! ^ R (^ M ^( R » 



dR / dR' e *CH R - R ') 



a -2ifc 8 [f(R)-£(R')] 



}-<< 



,-2jfc i {(R)^ e 2ifc^(R 



G 



Assuming translational invariance again and specializing to the specular direction, one obtains with a second order 
cumulant expansion (see Ref. J33[ for a more detailed derivation) 



P(0) 



-0E 



1 + i J dR - 1 



(19) 



where for brevity: 



= (2a) 2 , E = fc 2 . 



(It was not sufficient to simply set Q = in Eq.fllq) since then the integral would diverge due to the asymptotic 
properties of C(R).) 

I now wish to evaluate the specular intensity I for some specific correlation functions. To be consistent with the 
truncation of the cumulant expansion at second order, I will investigate the case of short-range order, in which higher 
order moments of the corrugation function do not play a significant role. One is thus led to consider two types of 
short-range correlation functions. These are of course just convenient models, assumed for lack of further knowledge 
of the probability density distributions. In both cases I will be concerned with evaluating the integral term in Eq.(|19[), 
denoted by: 



F(E) 



1 

27 



dR (e^W - 1 



1. Gaussian Correlation Function 

Let us assume a Gaussian form for the correlation function: 

C{R)=e-W>\ 

where I is the correlation length. Then due to cylindrical symmetry 



F(E) = 



l r 



RdR, 



where L is the linear extent of the surface. As such the integral cannot be evaluated analytically, but its derivative 
with respect to the energy can: 



OF _ I 2 
dE~ 2E 



e 0E _ e (3Ee- {L ' 1 ^ 



(20) 



Since short-range order was assumed one may safely neglect exp (— (L/Z) 2 ). Differentiating Eq.(|l9|) and combining 
with Eq.@ then yields: 

dE + P I + -—E— = °- 

This equation can be used conveniently for a best-fit of the correlation length I and the variance a 1 — /3/4, by using 
the experimental specular data, as a function of incidence energy E = fc 2 . 



2. Exponential Correlation Function 



Next assume a longer range, exponential form for the correlation function: 

C{R) = e^ R ' l \ 

where I is a new correlation length. Then: 



F(E) = 



,/3B(e-( R /"-l) 



1 



RdR. 
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Again the integral cannot be evaluated analytically, but its derivative can: 

BF I 2 

— = - [m((3E) - log(PE) - 7] . (21) 

[neglecting exp(—L/l)]. Here Ei(x) = e x (l/x + J °° e~ t (x — t)~ 2 dt) is the exponential integral and 7 is Euler's 
constant. Differentiation of Eq.([lJ) in combination with Eq.(|2l]) now yields: 

which can be used to best-fit the correlation length and the variance for a model of exponentially decaying correlations. 



IV. SUMMARY AND CONCLUSIONS 



This paper has presented new results on the inversion of structure of randomly corrugated surfaces from atom 
scattering data, within the SA. This work has been largely formal, with applications to be worked out in the future 
in connection with presently unavailable experimental data. The analysis presented here showed that within the 
framework of the SA, the scattering intensities in principle contain the full statistical information characterizing 
the surface disorder. Several potentially useful expressions were derived from which statistical parameters can be 
extracted in simple He scattering experiments. One application of the theory presented here is contained in our work 
on scattering from Ag/Pt(lll) |L3|, in which it was shown how randomly distributed adsorbates affect the scattering 
intensities. Additional theoretical applications will be undertaken in the future, but it is hoped above all that this work 
will stimulate experimentalists to further utilize inert atom scattering in the study of increasingly complex surface 
disorder. The results presented here suggest that such experiments can reveal a wealth of information concerning 
disordered surface structure, in particular on the statistics of randomly corrugated surfaces. 
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